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1.  Introduction  and  Summary 

In  this  paper  we  wish  to  establish  a  solution  procedure  designed  to 
solve  a  model  which  seeks  to  maximize  a  measure  of  asset  value  of  owned 
resources.  These  resources  are  updated  over  time,  and  change  according  to 
the  amount  of  manpower  and/or  maintenance  supplied. 

The  specific  model  which  we  address  is  a  slight  modification  of  a 
model  formulated  by  Rolf  Clark  in  a  pair  of  notes  ([1],  [2]).  It  is  a 
dynamic  model  with  a  system  of  difference  equations  linking  the  variables 
in  adjacent  time  periods.  The  objective  function  (to  be  maximized)  is  a 
discounted  sum  over  time  periods  of  a  measure  of  the  "effective  asset 
value"  of  the  owned  resources. 

Viewed  as  a  mathematical  optimization  problem,  the  model  appears 
to  be  a  difficult  nonconvex  optimization  problem.  In  Section  2,  we  will 
establish  the  notation  and  the  equations  describing  the  model. 
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The  nature  of  the  problem  is  such  that  we  were  able  to  accurately 
approximate  the  objective  function  and  design  a  branch  and  bound  method 
to  solve  the  problem.  The  details  of  this  "decoupling"  are  presented  in 
Section  3. 

A  moderately  expensive  way  to  implement  the  algorithm  suggested  by 
the  results  of  Section  3  would  be  to  design  (or,  perhaps,  modify)  a  code 
specifically  for  the  problem.  Fortunately,  we  were  able  to  effect  a 
solution  of  a  small  model  by  means  of  a  0-1  code  existing  at  GWU  (UNDO) . 
The  results  are  presented  in  Section  4. 

2.  The  Model 


The  model  presented  here  is  based  on  Clark Ts  notes  [1]  and  [2], 
and  also  on  a  Technical  Memorandum  by  one  of  the  author  [3]. 


State  Variables 


B 

t 

P 

t 

MXS 

t 

MBLOG 

t 


MX 


t 


MS 


t 


budget  at  time  t 
procurement  at  time  t 
maintenance  supplied  at  time  t 
maintenance  backlog  at  time  t 
maintenance  demand  during  year  t 
manpower  supplied  during  year  t 
manpower  demand  during  year  t 
asset  value  at  time  t 
assets  retired  at  time  t 
effective  asset  value  at  time  t 


The  state  variables  are  determined  by  a  system  of  difference 
equations, initial  values,  and  settings  of  the  control  variables. 
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Control  Variables 

Fl^  fraction  of  maintenance  supplied/maintenance  demanded 

F2^  fraction  of  manpower  supplied /manpower  demanded 

The  model  to  be  presented  has  a  number  of  parameters  which  are  set 
exogenously.  These  parameters  are  listed  below,  with  values  set  for  the 
illustrative  example  of  Section  4. 


Parameters 


c 

= 

1000 

al 

= 

0.03 

a2 

= 

2 

a3 

29/30 

a4 

= 

0.04 

a5 

= 

-0.01 

a6 

= 

25x10 

a7 

= 

0.5 

T 

= 

5 

\ 


► 


current  values 


y 

(number  of  time  periods) 


Initial  Values 

Bo  =  10 

AVq  =  100 
MBL0Gq  =  1 


The  relations  which  link  the  variables  over  the  time  periods  are 
now  presented.  The  basic  assumptions  underlying  and  defining  the  model 
are  exhibited  in  these  relations,  and  are  briefly  discussed  after  the 
relations  are  displayed.  The  reader  might  note  that  a  number  of  these 
relations  do  not  represent  constraints  on  the  problem  but  are  definitional. 
Indeed,  the  first  order  of  business  will  be  to  simplify  the  relations. 
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Relations 


The  maintenance  and  manpower  relations  are 


(1) 

B 

t 

= 

V1  +ai),: 

(2) 

P 

t 

= 

B  -  MXS  -  MS 
t  t  t 

(3) 

MXSt 

= 

Flt  (MBLOGt_1  +  MX  ) 

(4) 

MBLOG 

t 

= 

“2 

ao  MBLOG  ,  +  MX  -  MXS 

3  t-1  t  t 

(5) 

MX 

t 

= 

a. (AV  ,  -  a_  MBLOG  ,) 

4  t-1  7  t-1 

(6) 

MS 

t 

= 

F2  M 
t  t 

(7) 

M 

t 

= 

c(l  +  a5)t  (AVt_1)  a6 

(8) 

AV 

t 

= 

AVt-l  +  Pt  -  Rt 

(9) 

R 

t 

= 

AVJ;_1  (l-a3) 

(1) 

implies 

that 

the  budget  at  time  t  is  a  result  of  an 

initial  budget  growing  at  a  yearly  rate  of  .  Relation  (2)  defines 

procurement  at  time  t  to  be  the  difference  between  the  budget  at  t 
and  the  amounts  allocated  to  manpower  and  maintenance  during  time  period  t 

Relations  (3)  and  (6)  are  best  interpreted  as  defining  the  ratios 
FI  and  F2^_  as  the  fraction  of  maintenance  and  manpower  supplied  over 

the  amounts  of  these  quantities  demanded. 

The  maintenance  backlog  at  time  t  defined  by  (4)  is  interpreted 
as  a  fraction  of  the  backlog  in  the  previous  time  period  plus  the 
maintenance  demand  during  time  t  minus  the  maintenance  supplied  during 
this  time  period.  The  maintenance  demand  during  period  t  is  interpreted 
by  (5)  to  be  a  fixed  fraction  of  asset  value  minus  a  fraction 

of  the  maintenance  backlog  of  the  previous  time  period. 
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The  manpower  demand  as  defined  in  (7)  is  a  multiple  of  the  asset 
value  of  the  previous  time  period,  and  the  asset  value  in  (8)  is  updated 
from  the  asset  value  of  the  previous  time  period  by  adding  total  procure¬ 
ments  and  subtracting  assets  retired. 

Finally,  assets  retired  in  time  period  t  are  assumed  to  be  some 
fixed  fraction  of  the  asset  value  of  the  previous  time  period. 

The  reader  should  note  that  an  extension  of  this  model  which  allowed 
for  time-dependent  parameters  could  easily  be  incorporated,  if  meaningful 
parameters  could  be  estimated. 

We  are  now  in  a  position  to  simplify  the  relations,  by  eliminating 
variables  which' are  dependent  and  not  constrained. 

First,  we  may  let  MXS^_  and  MS^  denote  the  control  variables 
(instead  of  Fl^  and  F2^)  .  We  note  that  the  above  system  can  be 

simplified  considerably  if  we  let  MBLOG^  and  AV^  be  state  variables. 

For  we  may  rewrite  (4)  as 

a 

MBLOG  =  a„  MBLOG  .  +  a. (AV  -  a_  MBLOG  . ) 
t  3  t-1  4  t-1  7  t-1 

-  MXSt  by  (5) 

a 

=  (a3  -  a4a?)  MBLOG^  +  a4  AV^  -  MXSt 

and  we  can  rewrite  (8)  as 


AV  =  AV  n 
t  t-1 

+  (Bt  -  MXSt  -  MSt)  - 

(1-a3)AVt_i 

by 

(2)  and 
(9) 

=  a  AV 

!  +  Bn(l+ai)t  -  MXS 

-  MS 

+* 

by 

(1) 

In  order  to  simplify  notation,  we  set 

MXS  =  x 

t  t 

MS  =  s 

t  t 

AV  =  a 

t  t 

MBLOG  =  m 


-  5  - 
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and  get 


(10) 

a  = 
t 

a0a  -  x 

3  t-1  t 

-  st +  Bo(Uai> 

(11) 

m  = 
t 

a2 

Vt-l  +  fa3  ' 

■  VVVl  ‘  xt 

(12) 

0  £ 

Xt  i  a4St-l  + 

(l-a4o7)mt_1 

(13) 

0  i 

st  <  c(l-ha5)t 

a,,  a  , 

6  t-1 

where  (12)  and  (13)  come  from  the  restrictions  Fl^_,F2t  £  [0,1] 

together  with  the  defining  relations  (3),  (5),  (6),  and  (7).  We  can 
simplify  notation  further,  by  re-defining  the  coefficients  above  so  that 


(14) 

a 

t 

* 

kiVi 

-  X 

t 

-  st  +  B0(l+ai) 

(15) 

m 

_ 

kna  - 

+  k0m  - 

-  X 

t 

2  t-1 

3  t-1 

t 

(16) 

0 

< 

x  < 

k,  a  -  + 

k^m 

t  = 

4  t-1 

5  t-1 

(17) 

0 

< 

s  < 

t  — 

kr  a  -1 

6 1  t 

Note  that  (14)  and  (15)  define  the  state  variables  a^,...,a  and 

m^,...,m  as  linear  functions  of  the  control  variables.  Also  (16)  and 

(17)  represent  a  system  of  state  constrained  upper  bounds  on  the  control 
variables . 

Equations  (14)  and  (15),  inequalities  (16)  and  (17),  and  initial 
conditions  a^  =  AV^  ,  m^  =  MBLOG^  constitute  a  set  of  linear  constraints. 

We  are  now  in  a  position  to  address  the  objective  function. 

Empirical  evidence  has  suggested  that  a  meaningful  objective  might 


-  6  - 


have  the  form 
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T 

(18)  Z  w  EAV 
t=l  fc  fc 


where  are  preassigned  "weights"  and  EAV^  is  some  measure  of 

the  "effective  asset  value"  at  time  t  .  It  has  been  suggested  (Clark  [2]) 
that  EAV^  should  be  the  asset  value  AV^_  at  time  t  ,  diminished  by  some 

function  of  (a)  the  ratio  of  the  maintenance  backlog  over  the  asset  value, 
and  (b)  the  ratio  of  the  manpower  supplied  over  the  asset  value.  We 
assume  that 

EAVt  =  AV  min  {f  (MBLOG  /AV  )  ,  g(MS  /Mt>} 

or,  in  the  simplified  notation, 

EAVt  =  at  min  {fO^/a^  ,  g (s fc/£ tat_1> } 

where  f  and  g  are  to  be  described,  and  SL  =  c(l-Kx  )t  a.  (see 

t  ob 

equation  (7)). 

At  this  point,  we  wish  to  emphasize  that  the  exact  form  of  f  and 
g  is  immaterial  to  the  algorithm  which  follows.  The  algorithm  applies 
to  any  given  function  f  and  g  . 

A  function  f  which  has  a  form  compatible  with  the  current  problem 
( f (0 )  =  1  when  there  is  no  backlog,  f(l)  =  0  when  the  asset  value  equals 
the  backlog)  is 


(19)  f (v)  =  (1+5 8v) (1-v)58 
and  the  function  g  which  we  use  here  is 

(20)  g(r)  =  r5 (6-5r) 

(following  Soland's  note  [5]).  These  functions  are  depicted  in  Figure  1. 


-  7  - 
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f(v) 


1 


1 


r 


0 


1  v 


0 


1 


Figure  1.  Example  Functions  f  and  g 


In  summary,  the  model  which  we  address  has  the  form 
T 

maximize  £  w^  min  {f(mt/at)  ,  g  (s^/i  } 

xt’st  -  0 

subject  to  constraint  (14),  (15),  (16),  (17)  . 

Note  that  this  problem  is  equivalent  to  the  problem  P : 


maximize  v 


subject  to 


T 

v  <.  £  t 

t-1 


(21)  tt  <  wt  at  f(mt/at) 


)! t=i 


T 


(22)  tfc  £  wfc  afc  g(st/S,ta 


t  t-l' • 


and  constraints  (14),  (15),  (16)  and  (17). 


-  8  - 
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We  note  that  this  is  a  nonconvex  optimization  problem  as  the 
functions  on  the  right-hand  side  of  the  above  inequalities  are  generally 
not  concave. 

3.  Solving  the  Model 

In  this  section,  we  define  an  approximation  to  problem  P  which  can 
be  solved  globally,  and  indicate  a  method  that  could  be  used  to  solve 
the  problem. 

We  first  address  the  function  f (v)  .  Let  0  =  v^  <  vj_  <  v2  < 

...  <  v  =1  be  a  set  of  n+1  points  partitioning  the  interval  [0,1] 


n 


into  n  subintervals.  We  wish  to  define  the  piecewise  linear 
approximation  to  f(v)  which  agrees  with  f  at  VQ’**',vn  9  anC^ 

linearly  interpolates  adjacent  values.  Figure  2  exhibits  the  approxima¬ 
tion  desired. 


f  (v)  ^ 


1 


Figure  2.  Approximation  of  f 


We  introduce  variables  0 


tO 


0 


tl 


,  • .  •  , 


0 


tn 


vn,v  ,  ...,v  ,  and  note  that 

u  l  n 


n 


n 


f(vt)  2 


-  9 


associated  with  the  points 


T-474 


where 


(23) 

V 

t 

n 

=  Z 
3=0 

0  .  v. 
tj  j 

(24) 

1 

n 

=  Z 

9t. 

- 

3=0 

3 

CD 

rt 

L-i. 

>  0 

(j-0,1, . .  .  ,n) 

and  at  most  two  ®tj»s  can  positive,  and  if  two  are  positive,  they 

must  correspond  to  adjacent  grid  points  bordering  v^  .  This  latter 
restriction  is  known  as  the  adjacent  weights  restriction  (AWR) ,  and  is 
used  in  [4]. 

We  wish  to  replace  the  function  f  in  problem  P  by  its  piecewise 
linear  approximation.  Thus,  the  term 

at f  (mt/  at^ 

becomes 

jo  l9tA]  f(V 


and  the  new  constraint  (23)  becomes 


m 


n 


—  =  Z  8  .v. 

at  j=0  J 


i.e.  , 


n 


m  =  Z  [0  .a  ]  v. 

jto  ^  J 


Setting 


(25)  z  .  =  0  .a 

tj  tj  t 


-  10  - 
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constraint  (21)  becomes 


(26) 

<  w 
—  t 

n 

Z  z 

(27) 

m 

n 

=  Z 

j=o 

z  .v. 

t 

j=o 

tj  3 

(28) 

a 

II 

z  . 

t 

j=o 

tj 

,  f (v . ) 
tj  y 


where  the  latter  constraint  results  by  summing  (25)  over  j  ,  and  using  (24) . 

We  use  similar  manipulations  to  handle  constraint  (22) ,  but  we 
need  to  make  an  additional  approximation.  We  replace  at_j_  by  a^_  in 

(22)  in  order  that  the  construct  parallel  to  (25)  will  result  in  linear 
equations  similar  to  (26),  (27),  and  (28). 

Let  0  =  rn  <  r.  <...<  r  =1  be  a  set  of  m+1  points 
U  1  m 

partitioning  the  interval  [0,1]  into  m  subintervals.  Introducing 
variables  corresponding  to  these  subdivision  points  for  each  time 

period,  we  have 


if 


g  (r  )  *  Z 

j=0 


to  .  r 


tjj  J 


r 

t 


Z 

j=0 


<0  .r. 
1 3  J 


1 


Z 

j=o 


<0  . 

tj 


>  0 

tj  = 


and 


to 


tO 


* 


oo  satisfy  the  AWRs. 
tm 


11  - 
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Then  constraint  (22)  becomes  (replacing  a^  ^  by  a^_) 


(29) 

t  < 

m 

w  Z 

w  . 

t  “ 

c  J-O 

tj 

(30) 

s  = 

m 

a  z 

w  . 

t 

c  j-o 

tj 

(31) 

a  = 

m 

Z 

w  . 

t 

j-o 

tj 

where  we  have  defined  w  .  =  a  a)  .  .  Note  that  the  w  .  satisfy  the 

tj  t  t]  1 3  J 

AWRs  if  and  only  if  the  cj  .  do. 

t  J 


Summarizing,  the  problem  which  we  seek  to  solve  is  Problem  AP 


maximize  v 

T 

subject  to  v  <  Z  t 
~t=l 


n 

(26) 

t 

£ 

w 

Z 

z  .f(v.) 

t 

t 

j-o 

rt 

(_i 

L_l 

n 

(27) 

m 

= 

Z 

z 

V. 

t 

j-i 

o  j 

n 

(28) 

a 

= 

Z 

Z 

t 

j-< 

D  ’ 

tj 

m 

(29) 

w 

t 

Z 

j-o 

m 

(30) 

s 

= 

SL 

Z 

w  .  r . 

t 

t 

j-o 

tj  j 

m 

(31) 

a 

= 

Z 

w  . 

t 

j-o 

tj 

(14) 

a 

t 

= 

ki 

at-l 

-  \  '  st  +  Bo(1'tai)t 

(15) 

m 

_ 

k0« 

a 

+  k„m  ,  -  x 

t 

2 

t-1 

3  t-1  t 

12  - 
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(16)  0  <  x  <  k.a  _  +  krm  _ 

=  t  =  4  t-1  5  t-1 


(1?)  o  <  st  <  k6tat_1 


over  nonnegative  variables  v,  t,z.,m,a,s,w.,  and  x  for 

t*  tj’  t*  t’  t’  tj’  t 

t=l,...,T,  and  j=l,...,n  (for  the  z  . f  )  and  j=l,...,m  (for  the  w  . f  ) 

tj  s  tj  s 


Also,  the  sets  of  variables 


{z  }T  and  {w  , ...,w  }T 

to  tn  tO  tm  ^ 


must  satisfy  the  AWRs. 

It  is  this  latter  restriction  which  allows  us  to  classify  the  above 
problem  AP  as  a  nonconvex  problem.  An  algorithm  [4]  for  a  somewhat 
similar  problem  was  developed  at  GWU,  but  assumes  that  the  variables 
which  are  required  to  satisfy  the  AWRs  must  sum  to  1.  Undoubtedly,  the 
methodology  of  that  algorithm  could  be  extended  to  encompass  our  problem, 
but  would  require  some  delicate  reprogramming  of  the  working  code,  and, 
although  it  would  be  a  very  useful  addition  to  our  current  tools,  it  would 
also  take  several  months  to  complete. 

In  order  to  illustrate  the  solvability  of  problem  AP,  we  decided 
to  employ  LINDO  (see  [6]),  which  solves  problems  with  0-1  variables,  and 
adds  constraints  to  insure  that  the  AWRs  are  satisfied. 

We  note  that  equation  (28) : 


n 

a  =  Z  z  . 
c  j=0  » 

must  be  satisfied  by  a  set  of  variables  {z  ^,...,z  }  which  satisfy  the 

tO  tn 

AWRs.  For  this  set  of  variables,  we  introduce  0-1  variables  y  y 

tl  tn 


and  the  restrictions 


(32) 


200y  .  +  z  .  -  +  z 

tj  t,j-l  t,j 


>  a 
=  t 


(3  1 » •  •  •  j  n) 


-  13  - 
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(33)  Z  y  .  =  n-1 


j=l 


tj 


y  .  =  0,1 

1 3 


Equation  (33)  guarantees  that  exactly  one  variable  y  will  be  0  , 

i.e.,  exactly  one  of  the  constraints  (32)  will  be  binding.  Because  of 

constraint  (28),  the  z  ,  not  involved  in  the  binding  constraint  will 

tj  s 

be  zero. 


4 .  Solution  of  an  Illustrative  Example 

We  solved  problem  AP  with  the  0-1  option  of  LINDO.  The  following 
data  were  used  on  the  example  problem. 

The  parameter  values  as  given  in  Section  2. 


(w1, 

w2, 

w3, 

V 

w5)  = 

(i, 

2, 1,1, 2) 

o 

1  > 

V 

V 

V 

V  - 

(0, 

0.0005,  0.05,  0.1,  1.0) 

O 

1  u 

V 

V 

V 

V  = 

(0, 

0.1,  0.5,  0.9,  1.0) 

The  solution  obtained  is  exhibited  in  Table  1,  and  a  listing  of 
the  program  used  and  all  solution  values  is  copied  in  the  appendix. 

The  problem  had  102  constraints  and  136  variables,  20  of  which  were  0-1. 
This  is  a  relatively  large  problem  for  LINDO,  and  took  5529  iterations 
and  377  branches,  involving  10613  pivots. 
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Table  1 


Time  Period 


Solution  of  Example  Problem 
Main  t  enan  c  e  Manp  owe  r 


Supplied  (x^) 


Supplied  (st> 


Asset 


Value  (at) 


Maintenance 
Backlog  (m^) 


1 

4.86 

2.48 

100.00 

0.05 

2 

3.91 

2.45 

101.29 

0.14 

3 

4.05 

2.46 

102.61 

0.13 

4 

4.07 

2.46 

104.12 

0.16 

5 

4.17 

2.47 

105.75 

0.00 

optimal  value  *  711.19 
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PROBLEM  LISTING 
and 


COMPLETE  SOLUTION 
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Problem  Listing 

M AX  E  AV 

SUBJECT  TO 


2  ) 

X 1  +  i 

Ml  -  4.9144 

3  ) 

XI  <  = 

4.98 

4  ) 

0 . 96667 

A 1  +  A  2  +  X  2 

=  8.4872 

5) 

~ 

0.04  A 1 

-  0.91444  HI  x 

X 2  +  M2  ~ 

0 

6) 

~ 

0.0  4  A  1 

-  0 . 9 S  HI  +  X 2 

<  -  0 

7  ) 

- 

0 . 96667 

A 2  +  A3  +  X  3 

-  8.7418 

8  ) 

~ 

0.04  A2 

-  0.91444  M2  + 

X3  +  M3  » 

0 

9  ) 

0.04  A2 

-  0.98  M2  +  X3 

<  -  0 

1  0) 

- 

0.96667 

A3+  A  4  +  X  4 

=  9.0041 

1  1  ) 

~ 

0.04  A3 

-  0.91444  M3  + 

X  4  +  M  4  - 

0 

1  2  ) 

~ 

0.04  A3 

-  0.98  M3  x  X  4 

<  -  0 

1  3) 

- 

0.9666/ 

A  4  +  A  5  +  X  3 

=  9.2742 

1  4  ) 

~ 

0.04  A  4 

-  0.09144  M4  + 

M  5  +  X  5  = 

o 

1  5) 

“ 

0.04  A4 

-  0.98  M 4  +  X5 

<  =  0 

16) 

0 . 0005 

Z 1 1  +  0.05  Z1 2  + 

0.1  Z  1  3  - 

M  1 

+ 

Z  1  4  « 

17) 

0.0005 

Z21  +  0.05  Z22  + 

0.1  Z 23  - 

M2 

+ 

Z  2  4  - 

1  S  ) 

0.0005 

Z-3 1  +  0.0  5  Z  3  2  + 

0.1  233  - 

M3 

T 

23  4  - 

1  9  ) 

0 . 0005 

1. 4  1  +  0.03  Z  4  2  + 

0.1  Z  4  3  - 

H  4 

+ 

Z  4  4  :-- 

20  ) 

0.0005 

Z  5 1  x  0.0  5  Z  5  2  + 

0 .1  2  5  3  - 

M  5 

Z54 

21  ) 

21  0  + 

Z 1  1  +  212+  2 

1  3  -  A 1  + 

Z  i 

4  — 

0 

22  ) 

Z  2  0  + 

Z21  +  Z22  +  2 

23  •-  A  2  + 

12 

4  :=: 

0 

2  3  ) 

Z  3  0  + 

£31  +  Z32  +  Z 

3  3  -  A  3  + 

2  3 

■4  “• 

0 

24) 

Z40  + 

Z  4  1  r  Z  4  2  +  Z 

4  3  -  A  4  + 

Z  4 

A 

0 

2  ) 

250  t 

Z  5  1  +  Z  5  2  +  X 

3  3  ~  A  3  + 

Z  j 

i 

"0 

26) 

200  Y 1 1 

+  Z  1  0  +  Z  1  i  -- 

H  1  -  0 

27  ) 

200  Y  1  2 

.  +  Z  1  1  x  Z  1  2  - 

A 1  >  -  0 

23) 

200  Y 1 3 

+  2  1  2  +  Z  i  3  - 

A  1  --  0 

2  9  ) 

200  Y 1 4 

+  Z 1  3  -  A 1  + 

Z  1  4  >  -  0 

30) 

Y  1  1  + 

Y  1  2  +  Y  1  3  +  Y 

1  4  -=  3 

31  ) 

200  Y 2 1 

+  7. 2  0  +  Z  2  1  - 

A  2  >■••  0 

3  2  ) 

2  0  0  Y  2  2 

+  22 1  +  Z22  - 

A  2  >-  0 

3  3  ) 

200  Y 2 3 

+  Z  2  2  +  Z  2  3  - 

A  2  > ~  0 

3  4  ) 

200  Y24 

+  Z23  -  A 2  + 

Z24  >-  0 

3  3  ) 

Y  2 1  + 

Y  22  +  Y  2  3  +  Y 

24  =  3 

3  6  ) 

200  Y 3 1 

+  Z30  +  23 1  - 

A3  >~  0 

3  7  ) 

200  Y32 

+  23 1  +  Z  3  2  - 

H  O  “  0 

38) 

200  Y  3  3 

+  Z32  +  Z 3 3  - 

A  3  >  -  0 

39) 

200  Y  3  4 

+  Z  3  3  -  A  3  + 

Z3  4  >  =  0 

4  0  ) 

Y  3  1  + 

Y  3  2  +  Y  3  3  +  Y  3  4  =  3 

41  > 

200  Y 4 1 

+  Z 4 0  +  Z41  - 

A  4  >  =  0 

42  ) 

200  Y 4 2 

+  24 1  +  Z 4 2  - 

A  4  >  -•  0 

43  ) 

20  0  Y 4  3 

+  24 2  +  Z  4  3  - 

A  4  >-  0 

4  4  ) 

200  Y  4  4 

+  Z  4  3  -  A  4  + 

Z  4  4  >  =  0 

45  ) 

Y  4  1  + 

Y  4  2  +  Y 4  3  +  Y 

4  4  --  3 

4  6) 

200  Y 5 1 

+  Z50  +  25 i  - 

A  5  >"  0 

4  7  ) 

200  Y  52 

+  Z  5  i  +  2  5  2  - 

A  5  y  —  0 

48) 

200  Y53 

+  Z  5  2  +  253  - 

A  5  >~  0 

49  ) 

200  Y  5  4 

+  Z  5  3  -  A  5  + 

Z54  >=  0 

50  ) 

V  3  1  + 

Y  5  2  +  Y  5  3  +  Y 

5  4=  3 

-  16  - 


0 

o 

0 

0 

0 
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Problem  Listing  -  Continued 


si ) 

LAV  -  T1  - 

1  iL  “* 

T  3  -  T  4  -  T  5 

<  =  0 

52  ) 

1 1  0  +  0.99953 

Z1  1 

+  0.19908  Z 1 2  + 

0.01509  Z 1 3  - 

( i 

53  ) 

2  Z20  +  1.99916  Z2 

1  +  0.39816  Z22 

+  0.03018  Z23 

-  T  2 

51  ) 

Z  3  0  +  0.99953 

Z  3  1 

+  0.1  9908  Z 3 2  + 

0.01509  Z  3  3  - 

T  3 

55  ) 

Z  4  0  +  0.99958 

Z  4  1 

+  0.1  9903  Z-4  2  + 

0.01  509  Z  4 3  - 

T  4 

56  ) 

2  Z50  +  1.99916  Z51  +  0.39816  Z52 

+  0.03018  Z  5  3 

-  T5 

5  7  ) 

A 1  +  X 1  = 

1  0  4 

.  9  1 

58  ) 

- 

A  1  +  U 1  0  + 

U1  1 

+  U  1  2  +  U  1  3  + 

U  1  4  =  0 

5  9) 

- 

A  2  +  U  20  + 

u  2 1 

+  1422  +  U  2  3  + 

U  24  =  0 

60) 

A3  +  U  3  0  + 

U31 

+  U32  +  U33  + 

U  3  4  =  0 

61  ) 

- 

A  4  +  U  4  0  + 

U  4  1 

+  U42  +  U  4  3  + 

U  4  4  =  0 

62  ) 

- 

A 5  +  U  50  + 

U  5  1 

+  U  5  2  +  U  5  3  + 

U5  4  =  0 

6  3  ) 

- 

T 1  +  0.00006 

141  1 

+  0.1  0938  141  2  + 

0.88  57  4  14 1  3  + 

U  1  4 

6  4  > 

T2  +  0.00011 

LJ  2  1 

+  0.21875  U22  + 

1.77147  U 2 3  + 

2  U2  4 

65  ) 

- 

T  3  +  0.00006 

U  3  1 

+  0.1  0938  (432  + 

0.88574  14  3  3  + 

U  3  4 

66) 

- 

T  4  +  0.00006 

144  1 

+  0.1  0938  (4  4  2  + 

0.88574  U  4 3  + 

U  4  4 

6?) 

- 

15  +  0.00011 

(45 1 

+  0.21  875  (452  + 

1.771  47  1453  + 

2  U  5  4 

68  ) 

0.1  u 1 1  +  0.5 

U1  2 

+  0.9  14 1  3  +  U 1  4 

40.404  51 

0 

6  9  ) 

0.1  U  2 1  +  0.5 

(422 

+  0.9  W  2  3  +  (4  2  4 

-  40.81219  S  2 

= 

7  0  ) 

0.1  U  3 1  +  0.5 

U32 

+  0.9  (433  +  (434 

-  41.2244  33 

= 

71  ) 

0.1  144 1  +  0.5 

U  4  2 

+  0.9  U 4 3  +  (4 4 4 

-  41.6408  84 

- 

72) 

0.1  U51  +  0.5 

(452 

+  0.9  U53  +  (454 

-  42.0614  S5 

= 

7  3  ) 

SI  <=  '  2.475 

74) 

- 

0.0245  A 1  +  S 

12  <- 

0 

75) 

- 

0.02426  A2  + 

S  3 

>  0 

*76) 

- 

0.02401  A 3  t 

S4  < 

-  "0 

77) 

- 

0.02377  A 4  + 

3  5  ■■■ 

>  0 

78  ) 

- 

A  1  +  U  1  0  + 

U1  1 

+  20  0  YU  11  >  = 

0 

79  ) 

- 

A 1  t  U 1  1  + 

U1  2 

+  2  00  YU  12  >-■ 

0 

80  ) 

~ 

A  1  +  U 1  2  + 

U  i  3 

+  200  YU  13  >- 

0  ' 

81  ) 

- 

A  1  +  14 1  3  + 

U1  4 

+  200  YU  14  >= 

0 

82) 

Y  U  1  1  +  Y  U  1  2 

+  Y  U  13  +  Y  U  1  4  - 

3 

83  ) 

- 

A2  +  1420  + 

U  2  1 

-r  200  YU21 

0 

84  ) 

- 

A2  +  142 1  + 

U22 

+  200  YU 2-2  >- 

0 

85) 

- 

A 2  r  U22  T 

U  23 

+  200  Y U 2 3  >  = 

0 

86) 

A2  +  1423  + 

U24 

+  200  Y U 2 >- 

0 

87) 

Y 14  2  1  +  Y14  22 

+  YU  2  3-  +  YU  2  4  = 

3 

88  ) 

- 

A3  +  1430  + 

U  3  1 

+  200  Y 14  3  1  >  = 

0 

89) 

- 

A3  +  14  3 1  + 

U  3  2 

+  200  Y U 3 2  >  = 

0 

90) 

~ 

A3  +  14  32  + 

U3  3 

+  200  YU33  >~ 

0 

91  ) 

- 

A4  +  1433  + 

'4  3  4 

+  200  Y  U  3  3  +  200 

Y 143  4  >=  0 

92  ) 

Y 14  3  1  +  Y 1432 

+  Y 

U  3  3  +  Y  U  3  4  = 

3 

93) 

~ 

A  4  +  14  4  0  + 

U4  1 

+  20  0  Y U  4 1  y  — 

0 

9  4  ) 

A  4  +  144  1  + 

U4  2 

+  2o0  YU42  >= 

0 

95) 

- 

A  4  +  1442  + 

U43 

+  200  Y  U  4  3  > - 

0 

96) 

- 

A  4  +  14  4  3  + 

(44  4 

+  200  YU 4 4  >= 

0 

9  7  ) 

YU  4  1  +  Y 14  4  2 

+  YU 43  +  Y U 4 4  - 

3 

98  ) 

- 

A  5  +  U  5  0  + 

1451 

+  200  Y  U 5 1  >  = 

0 

99  ) 

A5  +  145 1  f 

U5  2 

+  200  Y U 5 2  >  = 

0 

10  0) 

A5  +  U52  + 

U53 

+  200  YU53  > :— 

0 

101  ) 

- 

A  5  +  1453  + 

U5  4 

+  200  Y U  5  4  > - 

0 

1  02  ) 

YUS  1  +  YU 52 

+  Y  U  5  3  +  YUS  4  - 

0 

0 


•j 

0 


.  H  D 

N U  M BEK 


INTEGER  VARIABLES: 


20 


17 
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The  Complete  Solution 


OBJECTIVE  FUNCTION  VALUE 


1 ) 

71  i  .  1 93726 

VARIABLE 

VALUE 

REDUCED  LOST 

Y  1  1 

0 . OOOOOO 

-426.513763 

Y  1  2 

1 . OOOOOO 

0 . OOOOOO 

Y  1  3 

1  - OOOOOO 

0.0000  0  0 

Y  1  A 

1  -  0  0  0  0  0  0 

0.0000  0  0 

—  V  2 1 

1  -OOOOOO 

0 . OOOOOO 

Y  2  2 

0.000000 

-137. 946239 

Y23 

1 . OOOOOO 

0.000 0 0  0 

Y  2  4 

1 .0000 0  0 

0.000000 

Y31 

1 . OOOOOO 

0.000000 

Y32 

0.000000 

-204 . 508804 

Y  3  3 

1.000000 

0  .  OOOOOO 

Y34 

1  .  OOOOOO 

0 . 0  0  0  0  0  0 

Y  4  l 

1  .  0  0  0  0  0  0 

0.0000  0 0 

Y  4  2 

0 . OOOOOO 

•334.489441 

Y  4  3 

1 -OOOOOO 

0.000  0  0  0 

Y  4  4 

1 .OOOOOO 

0 . 0  0  0 00  0 

Y  *5  1 

0 . OOOOOO 

0 . 0  0  0  0  0  0 

Y52 

1  .  0  0  0  0  0  0 

0  -  0  0  0  0  0  u 

Y  5  3 

1  .  0  0  0  0  0  0 

0  -  00  00 0  0 

YU  4 

1  .0000  0 0 

0  -  0  0  0  0  0  0 

cm  i 

1 . OOOOOO 

,  ~  2  fe  ^  '0  2  V  V  V  j 

Y  U 1  2 

1  .000  0  0 0 

"/.O.  J.iT7  7  -J 

Y  U  1  3 

1  .  0  0  0  0  0  0 

»•••  r  ^  n  ^  n  r 

i  O  »  J  7  7  7  J 

Y  W  1  4 

0.00  0  0  0  0 

0  .  OOOOOO 

YU 2  i 

1 . OOOOOO 

-35.762081 

YU  22 

1 .OOOOOO 

-33 .762081 

YU  23 

1 .OOOOOO 

-55.762061 

CU24 

0.000000 

0.000000 

Y  U  3  1 

1 . OOOOOO 

0 .OOOOOO 

Y  U  3  2 

1 . OOOOOO 

0 . OOOOOO 

i  y  j  o 

1 . OOOOOO 

0.000000 

Y  U  3  4 

0.000000 

0 . OOOOOO 

Y  U  4  1 

1  .00000  0 

-24.913073 

Y  U  4  2 

1 .OOOOOO 

-24 . 9130  73 

Y  U  4  3 

1  .  0  0  0  0  0  0 

-24 .91 3073 

Y  y  4  4 

0.00 0 0 0  0 

0.000000 

y  u  3  i 

1 . OOOOOO 

- 5  7 . 05  97  65 

YU  5  2 

1 .OOOOOO 

-57.059765 

Y  U  ’j  3 

1.000000 

-57.059765 

YU  3  4 

0 . OOOOOO 

0.000000 

•  EAV 

71 1 . 1 93726 

0  .OOOOOO 
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The  Complete  Solution  -  Continued 


Z  1  0 

0 . 000000 

0.001067 

Z1  1 

100.045631 

0 . 000000 

Z1  2 

0 . 0  00  0 0  0 

2.026954 

Z 1  3 

0.00 0 0  0 0 

1 . 920272 

Z20 

0 . 000000 

0 .690080 

Z21 

99 . 568626 

0 . oooooo 

Z2  2 

1 . 71755 7 

0 . oooooo 

£23 

0 . 000000 

0.661322 

£30 

0.000000 

1 . 0230 6 1 

Z  3  1 

1  0  1  .  0  1  i  0  3  2 

0.000000 

Z32 

1 .595954 

0.00000 0 

Z  3  3 

0 . 000000 

0 . 980426 

Z  4  0 

0 . 000000 

1  .92341 9 

Z  4  1 

102.013145 

0  .  oooooo 

Z42 

- 

2.111281 

O',  oooooo 

Z43 

0 . 000000 

1  .843263 

Z50 

1 05.748553 

0.00000 0 

Z51 

0 . 000000 

0.0001 4 3 

152 

0.000000 

0.01 4265 

Z53 

0. 000000 

0.0285  '3  0 

Hi 

1 00 . 045631 

0.00  0  0 0 0 

XI 

4.864372 

0 . OOOOOO 

H  1 

0.050026 

0 . oooooo 

H  2 

101.286194 

0 . 0  0  0  0  0  0 

X  2 

3 . 9 1 2 1 10 

0.00 0 0 0 0 

i'i  2 

0 . 1  3  5  4  6  0 

0.00 0 o  0  0 

*.  -y 

H  O 

1 02 . 607 i 09 

0  .  u  0  0  0  0  0 

X  3 

4 . 045224 

0 . 0  0  0  0  0  0 

M3 

0 . 130316 

0 . 0  0  0  0  0  0 

A  A 

1 0 4.124428 

0  .oooooo 

X 

A . 066376 

0  .oooooo 

n  4 

0 . 1  5  n  5  7  I 

0 . oooooo 

f  i  •  j 

1  0  5 . 7  4  8  5  5  0 

0 . 0  0  0  0  0  «> 

i'i  'J 

0 . 000001 

0  .  0  0  0  0  0  0 

V  r 
/» 

4  .  1  7  ?  2  9  7 

0  m  0  0  0  0  0  0 

Z  1  4 

0.000 0 0 0 

0.00  0 v 0  0 

Z  24 

0 .000000 

0 . oooooo 

Z  3  4 

0 . 0  0  0  0  0  8 

0 .oooooo 

Z  A  4 

0 . 0000  0  0 

0 . oooooo 

Z  5  4 

0 . 000000 

0 .285299 

T  1 

99 . 993393 

0 . oooooo 

T  2 

199. 737488 

0 . oooooo 

7  3 

101. 286430 

0 . oooooo 

r  4 

1 02 . 39061 7 

0 .oooooo 

7  5 

20/  .  7 8 5 o  i 

0 . 0  0  0  'J  v  0 

w  1  0 

0 . 0  0  0  0  0  0 

0.00 0 0 0  0 

U1  1 

0.00000 0 

0.11 4210 

U  1  2 

0.0000  0 0 

0.461  9 5  0 

1 J  1  3 

0 . 4 57153 

0.000000 

W  1  4 

99 . 583478 

0 . oooooo 

19 


T-474 


The  Complete  Solution  -  Concluded 
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